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ABSTRACT 

The rate of structure formation in the Universe is different in homogeneous and clustered 
dark energy models. The degree of dark energy clustering depends on the magnitude of its 
effective sound speed and for Ceff = 0 dark energy clusters in a similar fashion to dark 
matter while for Ceff = 1 it stays (approximately) homogeneous. In this paper we consider 
two distinct equations of state for the dark energy component, Wd = const and Wd = wo + 

wi with Ceff as a free parameter and we try to constrain the dark energy effective 

sound speed using current available data including Snia, Baryon Acoustic Oscillation, CMB 
shift parameter {Planck and WMAP), Hubble parameter. Big Bang Nucleosynthesis and the 
growth rate of structures fas{z). At first we derive the most general form of the equations 
governing dark matter and dark energy clustering under the assumption that CeS = const. 
Finally, performing an overall likelihood analysis we find that the likelihood function peaks at 
Ceff = 0, however the dark energy sound speed is degenerate with respect to the cosmological 
parameters, namely 0^ and Wd- 
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1 INTRODUCTION 


We are living in a special epoch of the cosmic history where the 
expansion of the Universe is accelerated due to an unknown en¬ 
ergy component, usually dubbed dark energy (DE). This accelera¬ 
tion has been discovered observation ally using the luminosity dis- 
tance of Type la supernovae (Snlal llPerlm utter et ^ 1^^, 1 19981 . 
1 19991: iRiess et al.ll2^: lAstier et aljr2006l: Ijha et al.ll2007h . In ad¬ 
dition to this, other observations includ i ng the cosmic microwave 
background (CMB) jBenn etLet^alj2()03 ; Spergel et alj200i|. 200^; 

' ivi: 


IPlanck Collaboration XII]||2015l: Hanck_Coll^oratioi^2^^ 20U) 


[2004 Coleetal 

2005h and baryon acoustic oscillation (BAG) 

I Eisenstein et alj 

2 OO 5 I: ISeo & Eisensteinll2005l: iBlake et al.l 

2011b 


support an accelerated expansion. At a fundamental level there are 
two different approaches to describe the phenomenon of the cos¬ 
mic acceleration and indeed many efforts are devoted to investigate 
its deep nature both observationally and theoretically. One way is to 
consider a fluid with a sufficiently negative pressure dubbed DE and 
the other is based on the modification of the laws of gravity on large 
scales. The first approach comes in many different scenarios. The 
simplest one is a very tiny cosmological constant A in Einstein field 
equations that has a (negative) pressure equal to it s energy density 
and equation-of-state parameter Wd = — = — 1 ( IWeinber3 1 9891 : 
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ISahni & Starobins^ l200d : [Peebles & Ratral 1200 3l) . The overall 
theoretical cosmological model (cosmological constant plus cold 
dark matter to explain galaxy rotation curves and the potential 
well for structure formation) is called ACDM model. Despite being 
highly consistent with observational data, the ACDM model suffers 
of two theoretical problem s, namely the fine-tuning and the cos¬ 
mic coincidence problem jWeinberd 1 19891 : ISahni & StarobinsF] 
I 2 OOOI : [Peebles & Ratral [200 3h . Differently from the cosmologi¬ 
cal constant case with equation of state (EoS) u>d = —1, 

other dynamical models have been largely studied in the litera¬ 
ture and_usually_categorhed in_twobrancheSj_qimtessence mod¬ 
els iArmendariz-Picon^tj^ 200(F Cqpelandet an[2006h and k- 


models lArmendariz-Picon et al J[ 1 999[. [20001: [Cliiba et al.[ 

[2000L[2009[:lAmendola & Tsuiikawa[[2010ti . 

The simplest way to modify gravity is to consider Einstein- 
Hilbert Lagrangia n as a generic function of the Ricci scalar 
B (f(R) theories, [Schmidt 199^^Magnano&^okolow^ 1994 


Dobado & Marotd [l995[ : [Capozziello et al.[ [20031 : [Carroll et al.[ 

2004) or add extra-dimension m odels like in the DGP model 


I Dvali. Gabadadze & PorTati|[2000[) . Understanding which class of 
models is the real one is one of the biggest challenges for cosmol¬ 
ogy- 


In addition to the background evolution, large scale struc- 
tures provide valuable infor mation about the nature of dark energy 
jTegmark et alj[2004 . [2006h . Primordial matter perturbations grow 
throughout the cosmic history and their growth rate depends on the 
overall energy budget and on the properties of the cosmic fluids. 
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DE slows down the growth rate of large-scale structures. Struc¬ 
tures grow due to gravitational instability and DE acts opposing 
and reducing the growth rate. The growth rate of structures can be 
measured from the redshift space distortion (RSD). Inward pecu¬ 
liar velocities of large-scale structures generate a distortion that is 
directly related to the matter density contrast. 


ters, DE perturbations and Ce. Finally in section[5]we conclude and 
discuss our results. 


2 EFFECT OF DARK ENERGY SOUND SPEED ON THE 
GROWTH RATE OF MATTER PERTURBATIONS 


Since the cosmological constant does not change in space 
and time, it can not cluster like dark matter (DM) and it has a 
negligible contribution to the energy density of the universe at 
high redshift. On the other hand, dynamical DE can cluster and the 
amount of clustering depends strongly on its effective sound speed. 
The effective sound speed is defined as = Ce = (hereafter 
we use Ce) where Spd and Spd are the pressure and energy density 
perturbations for DE respectively and coincides with the actua l 
sound speed in the dark energy comoving rest frame (lHij|l99^ . 


In quintessence models we have Ce — 1 so 
not grow on sub-horizon scales while in 
effective sound speed can be tiny (Ce <C 1) 

DE perturbations can 
k-essence models the 

I Garriga & MukhanovI 

I 999 I: Armendariz-Picon et al.ll 19991. 12 OOOI; 

Babichev et alj 20061; 

Akhourv et alj|201lb and DE perturbations grow similarly to dark 
matter (DM) perturbations. The possibility of DE clustering has 

been studied by many authors (Erickson et al.ll2002l: Bean & Dord 

2 OO 4 I: Hu & ScrantonI 20041; Ballesteros & Riotto 20081; 

de Putter et alj 20 id; Sanone & Maierotto 

12012; iBatista & Pacd 

2 OI 3 I: Dossett & IshakI 20131; Basse et al 

2014i; Batistal 20l4 


Pace et alj|20l4 l Steigerwald et alj|2014h . In particular, it has been 


shown that the homogeneous DE scenario fails to reproduce the 
observed concentratio n parameter of the massive galaxy clusters 
jBasilakos et al.|[200^ . In this framework, Ide Putter et alj j2010l) 
pointed out that CMB and LSS slightly prefer dynami c al DE 
with Ce ^ I and r ecently iMehrabi. Malekiani & Pac3 l l2015tl 
and iBasilakoi l l2015tl have shown that clustering DE repro¬ 
duces the growth data better in the framework of the spherical 
collapse model. A simi lar conclusion was suggested also by 
iNesseris & Sapon^ l l2014ll . 

The growth rate / = is usu ally approximated by / = 

Dm as first introduced bv IPeebl^ ( Il993h . In this parametrization 7 
is the so called growth index and ca n be used to d i stinguish between 
DE a n d modified gravity model s llLinderll2004 iHut erer & LindeJ 
I 2 OO 7 I : iBasilakos & Pouri 12012l : iRapetti et alj 2013 ). It is well 
known that for a ACDM model 7 is independent of redshift and 
equal to 6/11. The evolution of the matter density Dm depends on 
the evolution of the Hubble parameter iT(a) and hence on the par¬ 
ticular cosmological model adopted. In this paper we consider two 
distinct models, a constant Wd and a dynamical Wd (a), and we con¬ 
sider Ceff as a free parameter. Then based on the linear regime we 
numerically solve the perturbed general relativity (GR) equations 
to evaluate the growth rate of matter in the presence of DE clus¬ 
tering. Using a Markov Chain Monte Carlo (MCMC) method we 
can constrain the cosmological parameters using Snia, BAO, CMB 
shift parameter, the Hubble parameter, the Big Bang Nucleosynthe¬ 
sis (BBN) and growth rate data fas{z). 


The structure of this paper is as follows. In section |2] we de¬ 
rive the equations governing the linear growth of matter perturba¬ 
tions in a general relativistic framework and show the effects of 
DE clustering on the growth rate of matter. In section [3 we present 
all the details of the observational data used in this work to con¬ 
strain the cosmological parameters including the DE sound speed 
and their uncertainties. In section |4l we provide for the first time 
(to our knowledge) an approximated solution of the growth index 
of matter fluctuations as a function of the cosmological parame- 


In this section we revise the fundamental equations necessary to 
our analysis. The sound horizon of DE with effective sound speed 
Ce in a FRW universe is given by: 



dx , 


( 1 ) 


where 'H = ^, the prime being the derivative with respect to con¬ 
formal time ( 77 ) and Ui an initial scale factor. The nominal Hubble 
parameter is given hy H — ^ and thus H — aH which implies 




( 2 ) 


where an overdot refers to a derivative with respect to the cosmic 
time (t). In the case of Ce — 1, pressure suppresses any DE per¬ 
turbation with the consequence that DE may cluster only on scales 
comparable to the horizon. 

The opposite situation holds if Ce 1. Indeed in this case 
DE can cluster in analogy to the DM component and perturbations 
will grow with time. DE clustering modifies the evolution of DM 
perturbation and thus it affects the rate of structure formation in the 
universe. 

We start our derivation of the relevant equations by consid¬ 
ering the line element of an expanding universe in the Newtonian 
gauge without anisotropic stress: 

ds^ = —(1 -I- 2(f>)dt^ + — 2(j>)dx^ , (3) 


where </ is the Bardeen potential. First-order Einstein equations in 
Fourier space are: 

+ {sn^ + k^) 4> = —^(D^<5m + Dd<5d)(4) 

</"-f = ^Dd^Jd, (5) 

\ a J 2 Opd 

where Dm = Ddm + Db (Dd = 1 — Dm) is the matter (dark 
energy) density parameter and 5m (5d) is the corresponding den¬ 
sity contrast. The first-order energy-momentum conservation equa- 
tions for a generic fluid with equation-of-state parameter w are 
jMa & Bertschinge^ll995t) 


S' 


e' 


-{i + w){e 


30') - 3 — 

a 



- — {l-2,w)e 


1 -I- ui 


e + 


-^k^S + k^<j>. 
l + w 


( 6 ) 

(7) 


These equations are correct for any fluid with p = wp (for dust 
w — 0 and for dark energy w = Wd), where 5 is the density con¬ 
trast, 6 is the div ergence of the fluid velocity (6 = ik'^Vi) and 
can be written as dSean & Dordl2004l) 

^ = c, + 3H{l + w){c,-cld)~ , (8) 

where = Ca is the DE adiabatic sound speed: 


Ca. = W - 


3H{l + w) ■ 


(9) 


Note that the second term in the right hand side of Eq. appears 
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because we demand pressure perturbations to be a gauge invariant 
quantity dSean & Dordl2004l) . For a perfect fluid, perturbations in 
the pressure are purely determined by the adiabatic sound speed 
but for an imperfect fluid dissipative processes generate entropic 
perturbations and therefore we have a more general relation. In 
this case, Ce acts like a proxy for pressure perturbations and the 
growth of perturbation in the DE component depends on the effec¬ 
tive sound speed and not on the adiabatic sound speed any more. 
In the following this statement will be confirmed by solving the 
perturbed equations numerically. 

To study the effect of the DE sound speed on structure 
formation, we consider a universe with pressure-less DM and 
a DE component with varying equation of state that we spe¬ 
cialise to tnd( 2 ) = Wo -f The latter parametrization is 

the well known Chevallier-Polarski-Linder (CPL) parametrization 
dChevallier & PolarskilbOoH : lLindeill2003h . We eliminate 6 from 
Eqs. lUl and Q and find two second order differential equations for 
the density contrast of DM and DE. In addition using ^ = otH-^ 

and = cd'tid these equations can be 
written in terms of the scale factor. Finally our desired equations 
governing the growth of DM and DE perturbations are: 


for the initial value of Sm and 


dSm , D r _ o 

, o “r -^m , H” — ‘-5m ^ 

cta^ da 

^+Ad^-fBd<5d = S'd, 

da^ da 

and the coefficients [see also Eq. are: 


a V ^2 


— 

5m = 0 


Sn. = 3^ -f - 

da^ a 




(11) 


( 12 ) 


2 + 




d(j) 


da a^'H? ’ 


Ad — 
Bd = 


W 


2 -f 3ca — 6wd 

/ '^1 

3 (Ce - Wd) ( 1 + ^ - 3Wd -1- 3Ca - 3Ce 


dwd 


5*d — (1 “b Wd) 

k^ 


■-b 


( W \d^ 


3 d(j) dwd 


1 -I- Wd da da 


where ^ (or -^) is a function of the scale factor and using Eried- 
mann equations we have 

n' 


1 Dm -b Dd(l + 3uid) 1 /-I I on \ 

w - -2 —— = - 5(1 


(13) 


These equation s are not in agreement with Eq. (44) in 
lAbramo et al.l ( l2009h . which were obtained in the limit of a mat¬ 
ter dominated universe = — i) and a constant Wd. To resolve 
this discrepancy, see appendix 

We integrate Eqs. l IlOI and dl lb numerically from Zi = 100 to 
z = 0, in order to obtain t he density con t rast o f DM and DE. We 
use the same procedure of lAbramo et alj (l2009h to find the initial 
conditions. In the matter dominated era (f)' ~ 0, so from Eq. @ we 
have; 


dm,i — —2<j)i ( 1 -b 


k^ 


3uy 


(14) 


dSn 


da 


2jd_ 


(15) 


for its derivative. For dd the init ial value is set usin; 
the adiabatic perturbations condition jKodama & Sasaka 1 198 
lAmendola & Tsuiikawall2oI^ . 


dd,i — (1 “b Wd)dm,i , 
and its derivative is set to 


dSd i /1 I \ 


dwd^ 
da ^ da 


(16) 


(17) 


According to the above argument, by fixing the initial condition of 
(j)i we have all the initial conditions. We set = —6 x 10“^ which 
corresponds to (5m = 0.1 at present time for k — 0.1hMpc~^. Our 
results are robust under small changes of the initial conditions, and 
we don’t worry about the exact values. (For (pi = —7 x 10“®, (5m 
reach to 0.01 at present time but fag differs less than 10 “'*%.) 

DE clustering affects the growth of matter perturbations 
through the change of the potential (p. As we noticed the amount 
of DE clustering is directly related to its effective sound speed. 
We restrict our analysis to the choice of k — 1/A = 0.1/iM pc~* 
(10) which corresponds to A = 10(i“*Mpc dZhang et al.ll2012h . Note 
that the power-spectrum normalization crs which is the rms mass 
fluctuation on a scale Rg = 8 /i“*Mpc corresponds to fc = 
0.125/iMpc“*. On the other hand it has been common practice 
to assume that the shape of the power spectrum recovered from 
galaxy surveys matches th e linear matter power spectrum shape on 
scales k ^ 0.15/iMpc“* dSmith et al.ll2()03l : iTegmark et al.ll2004l : 
IPercival et al.ll2007h . Obviously the choice of fc = 0.1/iMpc“* as¬ 
sures that we are in the linear regime. We find that small variations 
around this value do not really affect the qualitative evolution of the 
growth rate of clustering and thus of 7 (z)Q 

To compare these results with observations we calculate the 
growth factor f{z) = — the growth index 7 ( 2 ) = 

dn ^(z) using our numerical results. The growth index in the 
ACDM model is redshift-independent and approximately equal to 
7 = 0.55. To compare this model to observational data we need to 
evaluate f(z)ag(z) where crg(z) is the mass variance in a sphere 
of radius of 8 Mpc/h. The variance ag(z) can be written in terms of 
ag at present time as ag(z) = ag(z = Also, in order 

to treat ag = ag(z = 0) properly for the DE models we rescale 
the value of ag by ag = ag a- Regarding crs a we utilise 

as.A = 0.818 (0.30 /Dm)°’^^ provided by the Planck analysis of 


ISnergel et al.l 1 201_5|) and it is also in agreeme nt with the results of 
Planck 2015 dPlanck Collaboration XIIll2015tl . 

DE perturbations not only depend on the sound speed but also 
on the EoS Wd. In the limit uta —^ — 1 all DE perturbations are 
washed out due to the 1 -\- Wd factor in front of the source term 
in the evolution equation of (5d. To show how the DE sound speed 
affects the linear evolution of DM, we consider D„i = 0.28 and 
h = 0.7 in the wCDM model to evaluate 5d and Ad = the 
relative DE density contrast, for a few distinct values of the DE 
sound speed as a function of the EoS. In Eig. Q the density contrast 
of DE as a function of Wd at the present time is presented. The 
non-clustering case remains homogeneous but for small values of 
the DE sound speed, the density contrast grows while increasing 


‘ Since we are in the linear regime we verify that for different values of k 
the differences in fag are practically negligible (~ 10“®%). 
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k = 0.1hMpc'\w-|=0.0 



Figure 1. The density contrast of DE as a function of the EoS at the present 
time for four different values of the sound speed. The red solid cui've shows 
a fully clustering DE model with Ce = 0. The green dashed (blue dotted) 
curve is for Cg = 10“^ (cg = 10“"^). A non-clustering model with Ce = 1 
is shown by a black dashed-dotted line. 


k = O.lhMpc \w-|=0.0 



Figure 2. The relative density contrast of DE as a function of the EoS at 
the present time for four different values of the sound speed. Line style and 
colours are as in Fig. G). 

the EoS. In contrast to the non-clustering case, the fully clustering 
regime with Ce = 0 gives a maximum value for the DE density 
contrast. In Eig. (H the relative DE density contrast is shown as a 
function of EoS. The behaviour of this quantity is similar to that of 
the density contrast. 

As we stated the quantity faB,{z) is affected by DE clus¬ 
tering. To show how fas{z) changes with the DE sound speed, 
we evaluate Afas{z) = x 100 and A-fiz) = 

^ 100 as a function of the EoS parameter. In the previ¬ 
ous equations, h stands for homogeneous DE. For the growth rate, 
results at present time are presented in Fig. (0. As expected, the 
deviation increases by increasing the EoS and for Wd < —0.9 the 
difference is less than 1%. The relative difference between homo¬ 
geneous and clustering DE for the growth index Aj(z — 0) has 
been shown in Fig. 0. The difference between the homogeneous 
and the clustering DE models is also very small for Wd very close 
to the ACDM model. 


k = 0.1hMpc'\w.|=0.0 



Figure 3. The relative difference of fas at the present time as a function of 
EoS. Line style and colours are as in Fig. o. 


k = 0.1hMpc’\wi=0.0 



Figure 4. The relative difference of the growth index at present time as a 
function of EoS. Line style and colours ai'e as in Fig. ffl. 


3 OBSERVATIONAL CONSTRAINTS ON THE DARK 
ENERGY SOUND SPEED 

In this section we use current available observational data sets 
to constrain the cosmological background parameters and the DE 
sound speed. In this analysis we assume that the DE sound speed is 
constant in time, regardless of the particular equation-of-state pa¬ 
rameter adopted. Our cosmological model will be described by the 
following parameters: DmO (matter density), Duo (baryon density), 
h — Ho/WO (normalised Hubble constant), wo and wi (dark en¬ 
ergy equation-of-state parameters) and Ce (effective sound speed) to 
describe the dark energy perturbations. In our analysis we assume 
a flat universe so that Ddm -f Ob + Da = 1, hence the amount of 
dark energy is known from the knowledge of the matter and baryon 
density parameters. 

The first data set we consider is th e Snia distance module from 
Union 2.1 sample dSuzuki et al.ll2012h . This data set includes 580 
SnIa and its is given by: 

,2 ^ [Mth(2i) - Mob(Zi)]^ 

Xsn — 2 _^ - ^ - , (fo) 

i ^ 
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Table 1. The current available BAO data which we use in our analysis. 


z 

di 

Survey & References 

0.106 

0.336 

6dF ZBeutler et al. 20111 

0.35 

0.113 

SDSS-DR7 (Padmanabhan et al. 2012) 

0.57 

0.073 

SDSS-DR9 (Anderson et al. 2013) 

0.44 

0.0916 

WiadeZ (Blake et al. 20111 

0.6 

0.0726 

WiaaleZ (Blake et al. 20111 

0.73 

0.0592 

WiaaleZ (Blake et al. 20111 


where /rth(«) = 5 logn, [(1 + z) + po, /tq = 42.384 - 

5 logj^Q h and ct; are the corresponding uncertainties. Before finding 
the minimum of Xsn we can expand Xsn around /tq 

Xsn = ^ + 2BflO + CfJ^o , 


(19) 


where 
A = 

B = 

C = 


E [uth{u0 = 0 ) ~ 

I '• 

E luthiuo = 0) ~ Uoh] 

E 4 . 


Obviously, for /ro = —BjC Eq. i fT^ has a minimum, namely A — 
^. Now by defining Xsn = A — ^, we can use the minimum of 
Xsn which is independent of /tq in order to find the best values of 
the parameters. Of course both estim ators provide the same results 
dNesseris & PerivolaronoulosIbOOSh . 

The second data set we consider is the BAO sample which in¬ 
cludes 6 distinct measurements of the baryon acoustic scale. These 
6 data points and their references are summarised in Tab. ([TJ. 

To find the X rao we follow the same procedure as 
iHinshaw et al.l l l2013l) . So the Xbao is given by 

xIao = Y^CbIoY , 


where Y = (d(O.l) - di, ^ 

di, d(0.6) — ds, d(0.73) — de) and 

f s (^^drag) 


( 20 ) 

ir>rf(0.44)- 


d{z) = 


with 


Dv{z) 


.(«)= r 

Jo 


Csda 


( 21 ) 


( 22 ) 


a^H{a) ’ 

is the comoving sound horizon at the baryon drag epoch, Cs the 
baryon sound speed and Dv{z) is defined by: 


Dviz) = 


{l + z)^Dl{z) 


H{z) 


(23) 


and Da{z) is the a ngular diameter distanc e. We used the fitting 
formula for Zd from lEisenstein & Hul ( Il998h and the baryon sound 
speed is given by: 

Ca(a) = ^ = , (24) 

3(1 + 

= 2.469 X ( lHinshawetaIll2013h . 


The covariance ma trix Cg^Q in Eq. J20t was obtained by 
IHinshaw et Ii]( l2013ll 

/ 4444.4 0. 0. 0. 0. 0. \ 

0. 34.602 0. 0. 0. 0. 

0. 0. 20.6611 0. 0. 0. 

0. 0. 0. 24532.1 -25137.7 12099.1 

0. 0. 0. -25137.7 134598.4 -64783.9 

V 0. 0. 0. 12099.1 -64783.9 128837.6 

The position of the CMB acoustic peak provides a useful data 
to constrain dark energy models. The position of this peak is given 
by (Za, R,Zt), where R is the scale distance to recombination 

Da{z,) 


la = 


TT- 


TsiZt) 

R = ff^HoDAiz,), 


(25) 

(26) 

and rs(z) is the comoving sound horizo n defined in Eq. l|22t. In 
this case we used the formula for fromiHu & Sueivam^ ( 1 19961) . 
Eor the WMAP data set we have ( IHinshaw et alj201^ 


XcMB = 


and 


-'CMB 


Za- 302.40 \ 

R - 1.7264 
- 1090.88 / 

3.182 18.253 -1.429 

18.253 11887.879 -193.808 

-1.429 -193.808 4.556 


(27) 


(28) 


In addition to this data set the Planck data provide more accurate 
CMB data for which the position of the acoustic peak is given by 
dShafer & Hutereill2014ll 

/ Za- 301.65 \ 

XcMB = R- 1.7499 , (29) 

V - 1090.41 / 


and 


-'CMB 


42.7044 -418.36 

-418.36 57366.3 

-0.7820 -762.152 


-0.7820 

-762.152 

14.6995 


(30) 


(31) 


In both cases the Xcmb is given by : 

XcMB = XqmbC CMB Xcmb . 

A further data set used in this work is t he Hubble evolutio n 
data obtained from the evolution of galaxies dSimon et al.ll2005]) . 
We use the 12 available data points and the for this data set is: 

[H{Zi) - Hndff 


Xh - ■ 


(32) 


The Big Bang Nucleosynthesis (B BN) provides a data point 
dSerra et al.ll2009l : lBurles et alj|200lh which constrains mostly fig. 
The Xbbn is given by 

[nltd - 0 . 022)2 


Xbbn — 


(33) 


0.0022 

The final data set used is the growth rate data. These data were 
derived from redshift space distortions from galaxy surveys includ¬ 
ing PSCs, 2DF, VVDS, SDSS, 6dF, 2MASS, BOSS and WiggleZ 
and the data with their references are shown in Tab. We solve 
Eqs. dlOb and dllb numerically to find f(z)as{z) and compute Xts 
with 

Ifo-siZi) - f(T8,ohf 


2 

Xfs — 


= E' 


where we set 


(34) 
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Table 2. The fcrs{z) data points including their references and surveys. 


z 


Reference 

0.02 

0.360 ±0.040 

Hudson & Turnbull (2013) 

0.067 

0.423 ± 0.055 

Beutler et al. QODl 

0.10 

0.37 ±0.13 

Feix et al. QOl^i 

0.17 

0.510 ±0.060 

Percival et al. (20041 

0.35 

0.440 ± 0.050 

Sons & Percival 120091: Teamark et al. 120061 

0.77 

0.490 ±0.180 

Guzzo et al. 120081: Sona & Percival 120091 

0.25 

0.351 ±0.058 

Samushia et ^^20^) 

0.37 

0.460 ± 0.038 

Samushia et a^20^) 

0.22 

0.420 ± 0.070 

Blake etal. (20111 

0.41 

0.450 ± 0.040 

Blake etal. (20111 

0.60 

0.430 ± 0.040 

Blake etal. (20111 

0.60 

0.433 ± 0.067 

Tojeiro^^^ (20^1 

0.78 

0.380 ± 0.040 

Blake etal. (20111 

0.57 

0.427 ± 0.066 

Reid et al. (20121 

0.30 

0.407 ± 0.055 

Tojeiro^^^ (201^1 

0.40 

0.419 ±0.041 

Tojeiro et al. (20^) 

0.50 

0.427 ± 0.043 

Tojeiro^^^ (201^1 

0.80 

0.47 ± 0.08 
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Table 3. The best value parameters and their l-a unceifainty for the wCDM 
model. 


Pai'ameters 

Best (WMAP) 

Best (Planck) 

h 

u.oyoo_Q QQ3Y 

U. /UO4t_0 0012 

^DM 

n oq 7'Q”I"0.0027 
U.zz/o_Q 0029 

U.^OU1_0 0010 

Dg 

o o 
o o 

o o 

6 d 

+ 1 

o 

o 

o 

n n/io9+0 0003 
U.U452_0 0002 

WO 

n ri.iQt2+0 0144 

-0.943610 0141 

u.yyfo_0 0053 

Ce 

0. 

0.001 

0-8 

0.837 

0.829 


The overall likelihood function is given by the product of the 
individual likelihoods: 

£tot = £sn X £bao X £cmb x £,h x £bbn x £fs , (35) 

and the total chi-square Xtot is given by: 

Xtot = Xsn + XbAO + XCMB + xlf + Xbbn + Xfs • (36) 

We calculate the total chi-square Xtot and find the best value 
of the parameters with an MCMC algorithm. The number of de¬ 
grees of freedom is — n^t — 1, where N = 616 and nm is 

the number of the fitted parameters. The results of this analysis for 
the wCDM, w(t)CDM and ACDM are summarized in Tabs. ([4j 
and IS respectively. 

To compare the DE models we have computed the corrected 
Akaike information criterion (AIC) jAkai y il974[Sugiu^ll978tl 
which, in our case, due to N/ nnt > 40, is given by: 

AIC = Xmira + 2nflt ■ (37) 

A smaller value of AIC indicates a better model-data fit. Of course 
it is well known that small differences in AIC are not necessarily 
significant and therefore, in order to assess the effectiveness of the 
different models in reproducing the data, we need to estimate the 


Table 4. The best value parameters and their l-cr uncertainty for the 
w(t)CDM model. 


Parameters 

Best (WMAP) 

Best (Planck) 

h 

p 

o 

o 

1 + 

o o 

b b 

o o 

W 

00 o 

0 7070+0 0012 

U. /U/U_0 0013 

^DM 

n 9904+0 0028 
u.zzo^_0 0027 

0 2361 +0 0012 
U.ZC>Oi_0 0011 

Dg 

0 0474+0 0005 

U.U^/^_0 0005 

n 0481 +0 0003 
U.U45i_0 0003 

U?o 

-1.0176lg;gg2| 

-0.95204lg;gg“ 

Wi 

0.3289lg;g3®g 

-0.18512tg;g?g| 

Ce 

0.002 

0. 

<^8 

0.840 

0.829 


Table 5. The best value parameters and their 1- 
model. 

■a uncertainty for the ACDM 

Parameters 

Best (WMAP) 

Best (Planck) 

h 

C4 --H 

o o 
o q 
d d 
+ 1 
00 

o 

CD 

u./UDy _0 0010 

^DM 

n 99 C 1 +0.0030 

u.zzoi_0 0029 

0 oocrn + O.OOlO 
o.zooy _0 0011 

Dg 

0 04150 +O OOO 6 
U.U4OO_0 0005 

n 0481 +0 0003 

U.U45i_0 0003 

0-8 

0.839 

0.829 

model pair difference AAIC= AIC^ — AICa;. The higher the value 
of 1 AAIC|, the higher the evidence against the model with a higher 
value of AIC. With a difference | AAIC| ^ 2 indicating a positive 


evidence and |AAIC| ^ 6 indicating a strong evidence, while a 
value IA AIC I ^ 2 indicates consistency among the two models. 
The results of our analysis are the following: 

(i) Using WMAP data: 

• For the wCDM model, Xmin = 586.53, nm = 5, so 
AIC=596.53 

• For the w(t)CDM model, Xmtn = 585.32, nnt = 6, so 
AIC=597.32 

• For the ACDM model, Xmin = 589.22, nnt = 3, so 
AIC=595.32 

(ii) Using Planck data: 

• For the wCDM model, Xmin = 595.76, nnt = 5, so 
AIC=605.76 

• For the w(t)CDM model, Xmin = 595.50, rint = 6, so 
AIC=607.50 

• For the ACDM model, Xmin = 595.79, nnt = 3, so 
AIC=601.79 

Concerning the best value of the dark energy sound speed we find 
that it tends to zero but the corresponding error bars remain quite 
large within Icr. In particular Ce lies in the range € [0., 1]. 

In order to investigate the range of validity for Ce, in Figs. © 
and we provide the la and 2a contours of our analysis. Note 
that in both plots the upper panels are for wCDM in which we 
present the confidence levels in the (ce,Dm) and (ce,w) planes, 
where Dm = D^m + Dp. In the bottom panels of Figs. ej and 
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Figure 5. The Icr and 2a contours of Qm{wCDM), u;{wCDM), 
u;o{w(t)CDM) and tt;i{w(t)CDM) versus DE sound speed using WMAP 
data. The Icr and 2a contours correspond to ~Xb ~ X^ ~Xb ~ 

6.16. The green (red) area correspond to la {2a) using only fas data and 
purple (blue) show Icj {2a) using all data set. 


Figure 7. The fa8{z) quantity (using Planck data), for the best values 
cosmological parameters for the wCDM (green dot-dashed curve) and 
w(t)CDM (red solid curve) models. The ACDM model is shown by the 
violet short-dashed curve. 



Figure 6. Same as Fig. but using the Planck shift parameter. 


© the contours for wo and wi in the CPL model are shown with 
respect to the DE sound speed. From this analysis it becomes clear 
that there is a strong degeneracy between Ce and (fln^jW) which 
implies that all values in the interval 0 ^ Ce ^ 1 are acceptable 
within the Icr uncertainty. As we stated, with data used in this pa¬ 
per the error bar of DE sound speed is quite large. 

In Figs. Q and we present the quantity /erg (2) for our 
best value parameters by considering the Planck and WMAP data 
for the wCDM, w(t)CDM and the ACDM models, respectively. We 
also show the observational data points. In addition to this quantity 
in Figs. © and not the growth index for the best values of the 
parameters have been shown. Note that using Planck CMB data our 
likelihood analysis indicates that all three models are very close to 
each others. 0 

Previous works in literature tried to put constraints on the 
dark energy ef f ective sound speed Ce using different kind of data. 
Ide Putter et alj ( l2010h combination of CMB temperature power 


^ See the results of for the Planck case. 



Figure 8. The fo-s{z) quantity (using WMAP data), for the best val¬ 
ues cosmological parameters for the wCDM (green dot-dashed curve) and 
w(t)CDM (red solid curve) models. The ACDM model is shown by the vi¬ 
olet short-dashed curve. 


spectrum data, their cross-correlation with several mass-density 
tracers and the SDSS ERG auto-correlation function. Supemovae 
data were used to break degenerac ies with background cosmolog¬ 
ical parameters. iHannestad )set of Supe rnova data, LSS and 
CMB power spectra. Finally Ixia et al.l l l2008h performed a similar 
analysis for a single perfect fluid and a two-field Quintom dark en¬ 
ergy model with w = —1 crossing by analysing CMB anisotropy 
data, LSS and SNIa observational data. In all these studies, using 
a similar approach to the one used in this work, the authors reach 
our same results. While previous and current data can constrain at 
a good level the current equation-of-state parameter of the dark en¬ 
ergy component, the quality of the observations is unfortunately 
still not sufficient enough to put any constraint on the dark energy 
effective sound speed. Note however that this is also due to the neg¬ 
ligible contribution of dark energy at early times on one side, and to 
the fact that curre n t obse rvations favour w ~ —1. As pointed out 
by Ide Putter et al.l (I 2 OIOI), if one consid ers the case of early dark 
energy models jporan & RobberslI^Obl) where the contribution of 
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Figure 9. The growth index (using Planck data), for the best values cosmo¬ 
logical parameters for the wCDM (green dot-dashed curve) and w(t)CDM 
(red solid curve) models. The ACDM model is shown by the violet short- 
dashed curve. 



Figure 10. The growth index (using WMAP data), for the best values 
cosmological parameters for the wCDM (green dot-dashed curve) and 
w(t)CDM (red solid curve) models. The ACDM model is shown by the 
violet short-dashed curve. 

dark energy at early times, i.e. CB, is not negligible, then more 
stringent limits can be set on Ce. 


4 GROWTH INDEX ANALYTIC SOLUTION 


In section we investigated the evolution of the growth index by 
solving numerically the syste m of Eqs. dlOb and jl lb . Here our 
aim is to extend the work of iBasilakoj 1 2015li in order to provide 
a general 7 ( 2 :) approximated solution which can be used in studies 
of structure formation. On sub-horizon scales, namely ^ 3> H 
(or ^ T-L^), Poisson equation (see appendix iBl takes the form 


-- 2 <P = [Sin 


L^m + Od^d(l + 3Ce)] . 


(38) 


Under the above conditions, Eq. Oil becomes 


dS H-^ da ^ 2 3ce)Od5d] . 


( 39 ) 


In this framework, for 5d = 0 , the latter equation reduces to the 
well known scale independent equation which is also valid for the 
concordance A cosmology. 

Concerning the equation of state parameter, it is well known 
that one can express it in terms of the Hubble parameter dSaini et al.l 
120001: iHuterer & Tumeill200lh 


Wd{a) 



2 „ din H 
3 ^ da 


(Q<) 


(40) 


or 


d\nH _ H 
“ da ^ IP 


3 

2 


-ind(a)Hd(a) , 


(41) 


where Hm(a) = 1 - f2d(a) = a^iP{a) = H{a)/Ho. 

Now, substituting Eq. j41b and / = dXnSya/d\na into Eq. l l39b we 
obtain the basic differential equation which governs the growth rate 
of clustering 


+ /^ + ^2 ~ “^ 2 ^ 3” 3ce)AdHd] , 

(42) 

where Ad(a) = 5d/5m. To this end, changing the variables in 
Eq. ( 142b from a{z) to redshift = —(1 + z )~^and utilising 
f{z) = we arrive to 


-(1-1- 2 ) 72 ln(Hm) -h Hm -h 3tndHd 



(43) 


where 72 = d'yjdz and 

V( 2 ) = l+^J^Ad(2)(l+3Ce). (44) 

On the other hand, the parametrization f{a) — 
dln^m/dlna ~ has a great impact in cosmolog¬ 

ical studies because it can be used in order to simplify the 
numerical calculations of Eq. HI- Obviously, a direct integration 
gives 


5m (a, 7 ) 



{nl{u) - 1 ) 


(45) 


where a{z) = 1/(1 -|- z) and Oi is the scale factor of the universe 
at which the matter component dominates the cosmic fluid (here 
we use Oi ~ 10“^ or Zi ~ 10). Hence, the linear growth fac¬ 
tor normalised to unity at the present epoch is D{a) — 
Therefore, in order to proceed with the analysis we need to some¬ 
how know the functional form of 7 ( 2 ). Erom the phenomenological 
point of view we may parametrize 7 ( 2 ) as follows 


7 ( 2 ) = 70 + 'yiy{z) ■ 


(46) 


This equation can be seen as a first order Taylor expansion around 

some cosmological quantity such as a{z) and z. _ _ 

Recently, it h as been found dBasilakosI I 2 OI 2 I : 
iBasilakos & Pouril 2012 L and references therein) that for those y{z) 
functions which satisfy the condition j/( 0 ) = 0 [or 7 ( 0 ) = 70 ], 
the parameter 71 is written as a function of 70 . For example, 
at the present epoch [2 = 0 , 7 z( 0 ) = 712 / 2 ( 0 ), Xq = A( 0 ), 
Wo = Wd(0)], Eq. ( 143b is written as 

»;L°o + 3wo(7o - Dfldo + I - 
yz ( 0 ) In Hmo 


(47) 
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- 0.2 - 0.1 0 0.1 0.2 


Figure 11. The pair (70 , 71 ) as a function of Ajq . The solid and the dashed 
lines correspond to the w(t)CDM and wCDM dark energy models, respec¬ 
tively. The homogeneous case Ajo = 0 is shown by the solid point. For the 
cosmological parameters, we use the values of Table 3 and 4 (third column). 

w here = dy/dz . Note that a similar equation has been found 
in ISasilakosI j2015l) in the case of Ce = Wd with Wd = const. 
As it is expected, for the homogeneous DE case (Aa = 0, 
X = 1), we verify that the above formula boils down to that 
of IPolarski & Gannouiil ( l2008h for y(z) = Within this frame¬ 
work , assuming y{z) = 1 — a{z) = llBallesteros & Riottd 
200^, we fully recover re s ults in literature jlshak & Dos settll 20091 : 
Bueno Belloso et ahlbOl ll : Idi Porto et alj|2012h . Notice that below 
we focus on y{z) = 1 — a{z) = with 3 / 2 ( 0 ) = 1. The fact 
that f2d(«) — 0 at 2 1 implies that the asymptotic value of the 
growth index 700 = 70 + 7 i is not really affected ^ the dark en¬ 
ergy clustering. Therefore, plugging 70 = 700 — 71 □ into Eq. ( 147 1 
we can obtain the constants 70,1 in terms of (flmo, tco, Aao, Ce). 

In Eig. ( fTTt we present (70,71) as a function of Aao- The 
curves are constructed using the parameters from Table 3 and 4 
(third column) and they correspond to w(t)CDM (solid) and wCDM 
(dashed) models. We observe that for Aao > 0 the growth index 
starts to deviate from that of the ACDM model, namely 70 < 0.55 
and 71 > 0. In the case of Aao < 0 the value of 70 is greater than 
that of the homogeneous case (70 > 0.55). In this context, con¬ 
cerning the value of 71 we find that it becomes negative. Of course 
for Aao = 0 the pair (70,71) reduces to that of the homogeneous 
case (see solid points in Fig. GD), as it should. 


5 CONCLUSIONS 

To summarize, we study the impact of dark energy clustering on 
the growth index of matter fluctuations. Initially we provide the 
most general form of the equations governing dark matter and dark 

® Regarding the asymptotic value of the growth index we use 700 ~ 
3(w — l)/(6w — 5) for th e wCDM model (see iLinder & CahnI |2007| : 
iNesseris & Perivolaropoul^ |2008l) and 700 ~ 0.55-I-0.05[1 + ut( 2 : = 1)] 
for the w(t)CDM model iLinderll^OO^ . 


energy clustering within the framework of Ce = const. Then us¬ 
ing the well known equation of state parameters, namely 'Uid(2) = 
wq -\- w\z/[1 + z), Wd[z) = const and the current cosmological 
data we place constrains on the cosmological parameters, includ¬ 
ing that of the effective sound speed Ce- Although the likelihood 
function peaks at Ce ~ 0, which indicates that the dark energy 
component clusters in analogy to the matter component Ce ~ 0, 
the corresponding error bars are quite large within 1 — a uncertain¬ 
ties which implies that Ce remains practically unconstrained. We 
also compared our findings with previous work reaching the same 
conclusion that at the moment the quality of cosmological data is 
not sufficient enough to put constraint on the dark energy effective 
sound speed. Future cosmological data, based for example on Eu¬ 
clid, are expected to improve even further the relevant constraints 
on Ce and thus the validity of clustered dark energy will be effec¬ 
tively tested. Finally, we have derived a new approximated solution 
of the growth index in terms of the cosmological parameters, dark 
energy perturbations and Ce- 
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APPENDIX A: PROOF OF EQ. (12) 

We start with Eqs. @ and ([7]l. The term appears in both equa¬ 
tions but it behaves very differently in these equations. In the first 
equation we have 

A '1 

-Sn^S = -mCeS - 9^(1 -b Wd){Ce - Ca)e , (Al) 

op fc 

and on sub-horizon scale we can neglect the latter term (fc^ H^), 
but in Eq. 0 we have 

CeS + 3'H{1 -l-uid)(ce — Ca)8 , (A2) 

dp 

where the latter term can not be neglected. Differentiating Eq. 
with respect to conformal time we have: 

S” + Wdd + {1 + Wd)9' + SH'cJ (A3) 

+ SHceS' - 3H'wdS - SHw'dS 

— SHwdS' = 3wd4>' + 3(1 + Wd)(j)' . 

Now from Eq. Q 

t r 

8' = -U{\-3wa)8- 6>+fc \ -b3T((ce-Ca)6)-bfeV , 
1 -b U)d 1 -b Wd 

(A4) 

and from Eq. nil 

, ,,, <5' 3Wce5 , 3Uwd& 

8 = 3<p - — -—- b —-- . (A5) 

1 -b Wd 1 + Wd 1 -b Wd 

Substituting Eqs. ( IA4t and ( IA5b into Eq. l lA3b . we have a second 
order equation governing the evolution of DE. Changing the inde¬ 
pendent variable to the scale factor, the coefficients in Eqs. (112b can 
be retrieved. On the other hand if we consider ^ = Ce and ignore 
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the second term in Eq. jA2t . we find 


Ad = - 


Bd = 


Ti' 

2 + -T—r + 3Ce — QWd 

3 (ce - Wd) (1 + ^ - 3wd) + - 3o-^ 


Sd = 


(1 + Wd) 

t 


3 

da^ 




a 

d(j) dwd 


2+^-3wd 


d(f) 

da 


0^14? 1 + Wd da da 

which coincide with the values in lAbramo et all ( l2009h for Wd = 
const and ^ = —i (matter dominated). We notice that for 
Wd = Ce = Ca. = 0 the coefficients for matter density contrast 
are recovered. 


APPENDIX B: POISSON EQUATION 

On sub-horizon scales, the basic equation describing the evolution 
of linear matter fluctuations is 

+ 2H (t)5m H- e<l> = 0 . (Bl) 

In this context the P oisson equation in the Fourier space is written 
as iLima et all 19971) 

(j) = — 47rGa^((5p -I- 35p) . (B2) 

where 5p = 5pm -f Spd and 5p — Spm + Spd- Now using Spm = 0, 
Spd = CeSpd, &pm = pmSm, &pd = Pd<5d, and inserting the above 
quantities into Eq.| |B 2I . we arrive to 

-2 47rG[pin5m + (1 + 3Ce)pd5d] , (B3) 

or 

-^0= 3Ce)nd5d] . (B4) 

Utilising the above equations it is easy to check that 

+ 2_f/(t)(5m = —^+ 3Ce)] . (B5) 

Obvio usly for Ce = Wd = con st, the latter equation reduces to 
that of lAbramo et aD 1 I 2 OO 9 I) and iMehrabi et all l l2015h . Changing 
the variables from f to a we finally obtain Eq. l l39t . 

This paper has been typeset from a TgX/ UTgX file prepared by the 
author. 
























